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SECTION  I:  GENERAL  PLASMA  THEORY  AND  SIMULATION 


ION  ACCELERATION  IN  THE  PLASMA  SOURCE  SHEATH, 
WITH  AN  INITIAL  ION  VELOCITY 

C.K.  Birdsall  and  S.  Kuhn^ 


In  an  earlier  noce,^  we  used  a  plasma  source  sheath  and  found  a  number  of  useful  results,  for  a 
monotonically  decreasing  potential,  such  as  the  ion  velocity  at  the  plasma  equalling  or  exceeding  the 
ion  sound  velocity  (hence,  obviating  the  need  for  a  Botun  pre-sheath  field  -  or  simply  placing  such  at 
the  source),  and  a  maximum  value  for  the  ion  flux  into  plasma. 


In  this  note,  we  will  still  approximate  the  ions  by  a  cold  beam  but  now  add  an  initial  ion  velocity, 
in  order  to  get  around  the  dilemma  of  infinite  ion  density  and  zero  ion  velocity  at  the  source.  See  Fig.  1 
for  the  model.  Assuming  that  most  electrons  are  reflected  to  the  source,  the  electron  density  is  still 
given  by  a  Boltzmann  factor,  with  coefficient  of  n„  ~  2n^ ,  where  the  -t-  means  that  injected  in  the  +  x 
direction.  The  ion  flux  density  is  given  by,  independent  of  dimension  z , 


r. 


=  Hi(x)  v.(z)  =  n* 


(1) 


We  take  the  ion  initial  velocity  to  be  given  by 


where 


(2) 


M 


Hence,  the  ion  density  is  given  by 

=  -y--  = 

v,(z)  'i  n  Vi(z)  y  n 


V 


(3) 


M 


The  problem  now  is  tc  solve  Poisson’s  equation,  in  one  dimension. 


Frequent  guest,  from  Institui  fur  Theosetische  Physik,  Univeniut  Iimsbiuck.  Innsbruck.  Austh*. 
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Figure  1 .  Mode)  for  planar  source  sheath,  with  warm  electrons  injecKd, 
along  with  cold  ions  injected  at  initial  velocity  . 
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(4) 


d^(x) 

dx^ 


(ni(x)  -  n.(x)) 
Co 


We  will  use  the  normalized  variables 


_  ►  X  "^io 

^  4  — ■  “  *  ''»=■ 


kT. 


^Dto  n*,  '  T, 

where  kT„K2n^  e^).  Thus,  Poisson’s  equation  now  is 


Tl'(^)= 


-a 


+  e 


rt(i) 


(5) 


(6) 


The  variables  ti  and  a  (“neuralization  parameter”,  “injection  ion  richness”)  are  those  of  Kuhn 

2  ^ 
(1981),  whereas  his  normalized  distance  was  =  i/Xo„,  with  ^D.o  =  Veo  kT„/(n*  e^).  Kuhn^  did 

the  complete  numerical  solution  for  the  case  Tq  =  1,  for  which  we  will  find  simple  analytic  results. 

(In  the  previous  note,^  the  initial  ion  velocity  was  zero,  just  as  it  was  in  the  derivation  of  the 
famous  Child’s  Law  of  1911;  the  physics  is  that  the  initial  ion  velocity  is  much  smaller  than  the 
directed  ion  velocity  in  most  places  elsewhere.  We  avoid  the  zero  here,  with  an  initial  ion  velocity 
equalling  the  average  velocity  of  a  half  Maxwellian  at  temperature  ■) 

Our  object  is  to  solve  Poisson’s  equation  so  as  to  obtain  the  potential  at  the  inflection  point 
(where  q  =  0),  which  we  take  to  be  the  plasma  potential,  ,  as  done  by  Kuhn  and  others.  This  solu¬ 
tion  is  from  (6),  as 

(qp<0)  (7) 

"  to 


This  expression  gives  us  the  normalized  ion  density  at  the  source  (i.e.,  the  neutralization  parameter  a) 
which  we  take  to  be  the  maximum  ion  density  in  this  system. 

We  now  can  plot  the  plasma  potential  q^  in  terms  of  the  injected  ion  richness,  cx,  as  done  by 
2 

Kuhn  (numerically)  in  his  Figure  2(a),  for  large  negative  collector  bias.  We  plot  our  results  on  top  of 
his  figure  and  find  excellent  agreement  for  the  plasma  potentia]  up  to  ion  richness  factor  a  =  2.  Our 
result  (7)  has  a  maximum  ion  richness  value  of  a  =  2.09  at  a  potential  of  q^  =  -0.18.  Some  numben  of 
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Figure  2.  Plasma  potential  Tip,  mean  ion  velocity  <v,>  and  slowest  ion  velocity  Vicuaff  plasma 

edge  as  a  function  of  ion  richness,  a;  numerical  results  taken  from  Kuhn^.  The  dashed 
lines  show  the  analytic  results  of  this  note,  really  overlaying  Kuhn’s  results,  up  to  a  =  2. 


interest  are; 


at  a  =  I.O  ,  71-  =  -1.58  ,  — ^  =  1.38,  with  v„  defined  below 

v„ 

Let  us  now  look  at  the  ion  speed  at  the  “plasma”,  defined  as  above.  The  ion  speed  in  the  plasma 
region  comes  irotn 


1  —2e  <I>  7  2  f^ui 


(9) 


Let  us  define  “sound  speed”  as 


M 


so  that 


V.  _  2  (t,-7ni,,) 
vi  ~  7t(x„+l) 


(10) 


For  Kuhn’s  figure,  the  ions  and  electrons  have  equal  temperatures  at  the  source,  x,  =  1,  so  that 
v„  =  •^2JcT„/M  and  the  ion  speed  in  the  plasma  region  is  found  to  be 


where  is  obtained  from  (7).  This  value  is  also  plotted  on  top  of  Kuhn’s  results  and  is  observed  to 
be  very  close  to  his  cutoff  velocity  (lowest  ion  velocity),  and  a  1  .tie  lower  than  his  average  velocity. 

The  point  is  that  in  our  cold-ion  model 


for 


a 


<  2. 


In  the  warm-ion  model,  the  inequality  <Vi>pi,TO  >  Vjomd.  holds  true  for  almost  all  values  of  a,  as  can 
be  seen  from  Fig.  2  Hence,  the  ions  are  at  least  at  sound  speed  upon  leaving  the  planar  source  sheath, 
obviating  the  need  for  further  Bohm  acceleration  in  the  (coUsionless)  plasma,  preceding  a  collector 
sheath. 
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SECTION  H:  PLASMA-WALL  PHYSICS,  THEORY  AND  SIMULATION 


A.  Collector  and  Source  Sheaths  of  a  Finite  Ion  Temperature  Plasma 
L.  A.  Schwager  and  C.  K.  Birdsall 

B.  Effects  of  Secondary  Electron  Emission  on  the  Collector  and  Source  Sheaths  of  a  Finite  Ion 
Temperature  Plasma 

L.  A.  Schwager 

C.  Effects  of  Ion  Reflection  on  the  Collector  and  Source  Sheaths  of  a  Finite  Ion  Temperature 
Plasma 

L.  A.  Schwager 

Three  titles  and  abstracts  follow.  These  will  be  ERL  reports,  soon. 

D.  Vortex  Dynamics  and  Transport  to  the  Wall  in  a  Crossed-Field  Plasma  Sheath 
Kim  Theilhaber  and  C.  K.  Birdsall 


An  interim  report  has  been  issued  (ERL  Memo  UCB/ERL  M87/18  10  April  1987)  in  this 
area,  with  title  and  abstract  following. 
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Collector  and  Source  Sheaths  of 
a  Finite  Ion  Temperature  Plasma 

L.  A.  Schwager  and  C.  K.  Lirdsall 

A  bstract 

The  region  between  a  Maxwellian  plasma  source  and  an  absorbing  surface  is 
modeled  with  an  electrostatic  particle  simulation  and  with  a  kinetic  plasma-sheath 
model.  In  the  kinetic  model,  Poisson’s  equation  and  Vlasov  equations  govern  the 
velocity  distribution  of  the  ions  and  electrons.  Our  numerical  and  theoretical  results 
for  collector  potential  and  plasma  transport  agree  with  the  bounded  model  of  Em- 
mert  ti  al.  but  differ  somewhat  from  those  using  traditional  Bohm  sheath  analysis. 
The  plasma  source  injects  equal  fluxes  of  half- Maxwellian  ions  and  electrons  with 
specified  mass  and  temperature  ratios  and  is  assumed  to  have  a  zero  electric  field. 
Representing  the  potential  change  within  a  distributed  full- Maxwellian  source  re¬ 
gion,  the  source  potential  drop  depends  primarily  on  temperature  ratio  and  evolves 
a  few  Debye  lengths  from  the  source  to  neutraUze  the  injected  plasma.  The  plasma 
flows  to  an  electrically  floating  collector  where  the  more  fanailiar  electron-repelling 
collector  sheath  appears.  Profiles  of  potential,  density,  drift  velocity,  temperature, 
kinetic  energv'  flux,  and  heat  flux  are  shown  from  simulation;  all  compare  very  well 
with  theory. 
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Effects  of  Secondary  Electron  Emission  on  the  Collector  and 
Source  Sheaths  of  a  Finite  Ion  Temperature  Plasma 

L.  A.  Schwager 

Abstract 

The  region  between  a  Maxwellian  plasma  source  and  an  absorbing  surface  which 
emits  cool  secondary  electrons  is  modeled  numerically  with  dynamic,  electrostatic 
particle  simulation  and  theoretically  with  a  static,  kinetic  plasma-sheath  model. 
The  coefficient  of  secondary  electron  emission  is  varied  up  to  and  beyond  critical 
emission  which  causes  electric  field  reversal  at  the  collector.  Results  from  these 
models  agree  very  well  over  this  wide  range  of  emission  coefficients.  Increasing  the 
secondary  emission  coefficient  reduces  the  collector  potential  which  increases  the 
total  energy  flux  to  the  collector  while  decreasing  the  ion  energy  deposited.  In  the 
simtilation,  some  heating  of  the  secondary  electron  stream  is  observed  to  gradually 
evolve  over  many  Debye  lengths,  possibly  becaiise  of  a  beam-plasma  interaction. 
This  heating  increases  potential  fluctuations  but  causes  only  small  deviations  from 
the  predictions  with  static  theory. 
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Effects  of  Ion  Reflection  on  the  Collector  and  Source 
Sheaths  of  a  Finite  Ion  Temperature  Plasma 

L.  A.  Schwager 

Abstract 

The  region  between  a  Maxwellian  plasma  source  eind  an  absorbing  siirface  which 
reflects  a  fraction  of  the  incident  ions  is  modeled  numerically  with  dynamic,  elec¬ 
trostatic  particle  simulation  and  theoretically  with  a  static,  kinetic  plasma-sheath 
model.  The  fraction  (  of  ions  reflected  is  varied  from  0  to  0.6  which  generally  in¬ 
creases  both  the  potential  drop  from  the  source  to  the  collector  and  the  energy 
transported  to  the  collector  surface.  Results  from  both  models  agree  well  when 
the  fraction  reflected  is  less  than  0.4  for  full  energy  transfer  to  reflected  ions.  With 
larger  fractions  and  with  slightly  less  than  full  reflected  energy,  simulations  show 
an  lon-ion  two-streaming  interaction  which  slightly  reduces  the  collector  potential 
drop  and  decreases  the  ion  energy  deposited  on  the  collector  surface  relative  to 
predictions  from  the  static  theory.  According  to  theory,  for  a  deuterium-tritium 
plasma,  a  collector  material  causing  the  reflected  ion  fraction  to  be  ^  =0.2  with  full 
reflected  energy  increases  the  magnitude  of  collector  potential  by  12%  and  the  ion 
energy  deposited  by  6%  over  those  predicted  when  (  =  0. 
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Vortex  Dynamics  and  Transport  to  the  Wall 
in  a  Crossed-Field  Plasma  Sheath  -  K.  Theil- 
haber  and  C.K.  Birdsall. 


Abstract 

Results  of  numerical  simulations  of  the  time-dependent  behavior 
of  a  transversely  magnetized  plasma-wall  sheath  are  presented.  These 
simulations  have  been  conducted  with  the  aim  of  modelling  plasma 
behavior  in  the  vicinity  of  the  limiters  and  walls  of  a  fusion  device. 
The  two-dimensional,  bounded  particle  simulation  code  “£82”  has 
been  used  as  a  tool  for  the  investigation  of  these  edge  effects,  in  an 
idealized  geometry  which  retains,  however,  the  essential  features  of 
the  physics  of  the  edge  plasma.  The  simulations  have  revealed  that 
the  bounded  plasma  is  subject  to  the  so-called  “Kelvin- Helmholtz” 
instability^  an  instability  maintained  by  the  non-uniform  electric  field 
which  is  induced  by  the  presence  of  the  material  walls.  This  instability 
is  seen  to  saturate  into  large  and  stable  vortices,  with  e0/Tj  ~  1,  which 
exist  in  the  vicinity  of  the  walls,  and  drift  parallel  to  their  surfaces. 
An  important  feature  of  these  structures  is  that  they  continuously 
convect  particles  to  the  walls,  at  an  “anomalous”  rate  much  greater 
than  that  induced  by  coUisional  diffusion,  a  feature  which  seems  tied 
to  the  mutual  interaction  of  the  vortices.  In  the  code  “ES2”,  volume 
ionization  of  neutrals  has  been  modelled  by  a  uniform  electron-ion  pair 
creation  in  the  simulation  region,  and  this  results  in  a  steady  state, 
in  which  the  linear  edge  instability,  the  nonlinear  fluid  dynamics  of 
the  vortices,  and  the  nonlinear  dynamics  of  the  particles  scattered 
by  the  vortices  all  balance  each  other.  This  steady-state  but  non- 
equilibrium  configuration,  which  is  a  first  model  of  the  edge  behavior 
induced  by  the  boundaries,  is  conceptually  analogous  to  Rayleigh- 
Benard  convection. 


^S.  Chandrasekhar,  Hydrodynamic  and  Hydromagnetic  Stability,  Oxford,  at  the  Claren¬ 
don  Press,  1961. 
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Professor  N.  Sato,  Tohoku  University,  Sendai,  Japan;  seminar,  'Basic  Experimenu  of  Plasma  Potential 
Formation  in  Magnetic  Minor  Formation,”  January  24,  1987 


Large  Electric  Fields  in  a  Magnetized  Plasma  Sheach;  Long-lived  Vortices  - 
K.  Theilhaber  and  C.  K.  Birdsall,  University  of  California,  Berkeley 

We  consider  Che  equilibrium  and  time-dependent  behavior  of  the  plasma-wall 
sheaCh  which  Is  formed  when  Che  external  magnetic  field  is  parallel  to  the  wall. 
The  presence  of  Che  material  boundary  initially  creates  a  layer  of  charge 
depletion  near  the  surface,  of  thickness  comparable  to  the  Ion  gyro-radius. 

This  results  in  a  strongly  nonuniform  electric  field  which  induces  a  sheared 
pCB  motion  of  Ions  and  electrons  parallel  to  the  wall.  Physical  arguments 
as  well  as  analytical  and  numerical  studies  suggest  chat  this  flow  is  subject 
to  the  Kelvln-Helmholtz  instability,  with  a  long-term  evolution  into  a  set  of 
vortices  separating  Che  boundary  fro  i  Che  plasma.  We  present  an  analytical 
model  for  Che  instability,  and  the  results  of  two-dimensional,  electrostatic 
particle  simulations  of  this  configuration,  in  Che  form  of  a  movie,  showing 
initialization  of  E,  Che  subsequent  shear  motion,  and  Che  buildup  near  Che 
wall  of  vortices  which  live  (at  least)  some  10s  of  ion  gyroperiods. 

DOE  Contract  FG03-  86ER53220  and  ONR  Contract  N00014-85-K-0809 
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Transport  induced  by  a  Crossed-Field  Sheath'  -  K.  THEILHABER, 
Plasma  Theory  and  Simxilation  Group,  University  of  California,  Berkeley, 
CA  94720. 


We  are  studying  transport  across  the  magnetic  field  in  a  plasma  bounded 
by  omducting  boirndziries  with  the  aim  of  modelling  plasma  behavior  in  the 
vicinity  of  the  limiters  and  walls  of  a  fusion  device.  Our  approach  is  to  use 
the  two-dimensional,  bounded  particle  simulation  code  ‘‘ES2’'  as  a  tool  for 
the  investigatioii  of  edge  effects,  in  an  idealized  geometry  which  retains, 
however,  the  essential  features  of  the  physics  of  the  edge  plasma.  Thus  far, 
our  simulations  have  revealed  that  the  bounded  plasma  is  subject  to  the 
so-called  “Kelvin- Helmholtz”  instability^  an  instability  maintained  by  the 
sheared  electric  field  which  is  induced  by  the  presence  of  the  material  walls. 
This  instability  is  seen  to  saturate  into  large  and  stable  vortices,  which  can 
be  described  by  some  variant  of  the  so-called  “Hasegawa-Mima”  equations®. 
These  vortices  exist  in  the  vicinity  of  the  walls,  and  drift  parallel  to  their 
surfaces.  An  important  feature  of  these  structures  is  that  they  continuously 
convect  particles  to  the  walls,  at  a  rate  much  greater  than  that  induced  by 
collisional  diffusion.  This  property  seems  tied  to  the  mutual  interaction  of 
the  vortices.  In  the  code  “ES2”,  we  have  modelled  volume  ionization  of 
neutrals  by  a  uniform  electron-ion  pair  creation  in  the  simulation  region, 
and  this  results  in  a  steady  state,  in  which  the  linear  edge  instability,  the 
nonlinear  fluid  dynamics  of  the  vortices,  and  the  nonlinear  dynamics  of 
the  particles  scattered  by  the  vortices  all  balance  each  other.  We  consider 
this  configtiration  a  first  crude  model  of  the  edge  behavior  induced  by  the 
boimdaries.  A  movie  will  be  shown,  shovnng  vortex  formation  and  the 
time-de|>endent  depletion  of  the  density  profiles. 
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'Thia  work  supported  by  U  S.  Deputment  of  Energy  Contract  No.  FG03-86ER53220. 
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don  Press,  1961. 
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PERFORMANCE  AND  OPTIMIZATION  OF  DIRECT  IMPLICIT  TIME 
INTEGRATION  SCHEMES  FOR  USE  IN  ELECTROSTATIC  PARTICLE 
SIMULATION  CODES* 

Richard  J.  Procasslnl  and  Charles  K.  Birdsall 
Electronics  Research  Laboratory 
University  of  California,  Berkeley 
Berkeley,  CA.  9‘i720 

and 

Bruce  I.  Cohen 

Lawrence  Livermore  National  Laboratory 
University  of  California 
Livermore,  CA.  9^550 


Implicit  time  integration  schemes  allow  for  the  use  of  larger  time 
steps  than  do  conventional  explicit  methods,  thereby  extending  the 
applicability  of  kinetic  particle  simulation  methods.  We  will  present  the 
results  of  a  study  of  the  performance  and  optimization  of  two  such  schemes 
which  are  used  to  follow  the  trajectories  of  charged  particles  in  an 
electrostatic,  particle-ln-cell  plasma  simulation  code,  the  two  direct 
implicit  schemes,  known  as  the  C^  and  D^  schemes,  were  chosen  for  their 
desirable  properties  which  allow  for  the  simulation  of  phenomena  at 
frequencies  much  lower  than  the  plasma  or  cyclotron  frequencies.  These 
properties  include  i)  the  relaxation  of  uAt  constraints  on  the  stability  of 
the  method,  li)  strong  damping  of  high-frequency  modes  for  which  uAt  >  1, 
and  lii)  second-order  accuracy  of  the  methods  in  the  simulation  of  low- 
frequency  phenomena  for  which  uAt  <  1. 

The  formulation  and  optimization  of  the  direct  implicit  schemes  will  be 
outlined,  followed  by  a  short  discussion  on  the  incorporation  of  these 
schemes  into  the  electrostatic  particle  code.  The  results  of  a  numerical 
stability  study  will  also  be  presented.  This  study  was  based  upon  the 
freely  expanding  plasma  slab  model,  which  was  run  for  various  combinations 
of  the  important  simulation  parameters  u  At,  Ax/i-  in  the  range  At 

i  200.0,  0.1  ^  ^  100.0,  excluding^fhe  region  for  which  v  At/Ax 

The  degree  of  numerlSal  heating  or  cooling  (lack  of  energy  consirvatlon)  was 
used  as  a  quantitative  measure  of  the  numerical  stability  of  the  scheme. 

The  results  indicate  that  the  D  scheme  possesses  a  much  larger  region  of 
parameter  apace  over  which  energy  is  well  conserved,  than  does  the  C. 
scheme.  A  locus  of  stable  operating  points  is  suggested  by  the  data  ror 
each  of  the  schemes. 


♦This  work  is  supported  by  the  USDOE  under  Contract  W-7M05-ENG-48. 
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Collector  Sheath  and  Source  Sheath  in  a  CoUicionleM  Finite  Ion  Temperature 
Plasma  with  Secondary  Electron  Emission  and  Ion  Reflection  at  the  Bounding  Sur¬ 
face  *  -  LOU  ANN  SCHWAGER,  PtoMtna  Theory  *»i  Simulaiion  Group,  Sleeironte*  JUeeurek  Let, 
Univereitp  of  Culifomia,  Berkeley  CA  94710. 


For  an  electrically  floating  surface  in  contact  with  a  plasma,  our  analytical  and  numerical  results 
for  the  surface  potential  and  plasma  transport  differ  considerably  from  those  usiag  traditional  Bohm 
sheath  analysis.  The  finite  sheath  region  between  a  coUisional  warm  plasma  and  the  confining, 
interacting  surface  is  modeled  numerically  with  electrostatic  particle  simulation  using  PDWl  and 
analytically  with  a  plasma-sheath  equation.  The  model  and  assumed  boundary  conditions  along 
with  a  typical  potential  profile  are  shown  below  For  the  purely  absorbing  surface,  we  derive  the 
plasma-sheath  equation  for  the  boundary  conditions  shown  below.  The  collisions]  plasma  source 
can  have  arbitrary  ion/electron  ratios  of  mass  and  t^nperature.  Our  technique  extends  the  analysis 
of  Kuhn  ^  with  a  temperature  ratio  Ti/T,  =  r  =  1  to  include  an  arbitrary  i«lue  for  r. 

Simulation  results  after  equilibration  and  our  time-independent,  analytic  theory  exactly  agree. 
When  compared  with  works  such  as  that  of  Chodura  ^  and  Stangeby  *  who  rely  on  a  Bohm  type 
of  semi-infinite  sheath  model,  our  derivation  predicts  a  lower  wall  potential  for  a  mass  ratio 
of  1836  and  r  =  1;  this  difference  is  within  20%.  Our  found  with  r  s  1  and  the  derived 
by  Stangeby  and  Chodura  W  r  =  0  agree  to  within  2%.  For  r  >  1  our  results  for  the  ion  drift 
velocity  vp  entering  the  collector  sheath  are  within  4%  of  Stangeby’s  assumption  that  this  velocity, 
normalized  by  the  ion  thermal  velocity  ot  is  (1  -t-  For  r  <  1  our  results  show  a  much  smaller 
vpi  e.g.  with  r  =  0.1  and  a  mass  ratio  of  1836,  «p/«t  s  1.8.  Note  that  the  presheath  region, 
which  we  define  below  by  4n ,  has  a  zero  electric  field;  hence  the  “preaheath  acceleration”  of  the  ioiut 
occurs  entirely  in  the  source  sheath  within  a  few  Debye  lengths  lengths  from  the  plasma  source.  Our 
potential  results  for  a  large  range  of  mass  and  temperature  ratios  compare  well  with  Emmert  '  who 
analytically  models  the  plsama-sheath  equation  for  a  bounded  system  with  Boltzmann  electrons. 

For  electron-induced  secondary  electron  emission,  we  have  simulated  emisaion  beyond  field  re¬ 
versal  at  the  plate  to  emission  cr^dents  F  over  100%.  For  emission  values  up  to  the  critical 
coeffident  at  field  reversal,  our  dependence  of^wonFforrsl  compare  well  to  that  of  Hobbs 
and  Wesson  '  who  assume  r  <  1 .  For  the  reflection  of  plasma  ions,  we  have  simulated  reflection 
coeffidents  A  up  to  80%.  The  wall  potential  decreases  as  we  increase  R.  Despite  an  ion-ion  two- 
streaming  interaction,  which  generates  substantial  potential  fluctuations,  the  time-averaged  values 
for  and  transport  agree  well  with  our  time-independent  analysis  which  imimn  r  1. 
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*  This  work  supported  by  U.S.  DOE  Contract  No.  DE-FG03-86ER53220. 

*  S.  Kuhn,  Plmsmu  Pkyhee  33,  886  (1981). 

*  R  Chodura,  yearns/  of  Nuclear  Material  111-112,  420  (1982). 

*  P.  C.  Stangeby,  Pkyt.  Fluids  27, 684  (1984). 

*  G.  A.  Enomert,  et  al,  Pkya.  Fluids  33,  806  (1980). 

*  G.  D.  Hobbs  and  J.  A.  Wesson,  Plasma  Physics  0,  85  (1967). 
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Large  Electric  Fields  In  a  Magnetized  Plasma  Sheath; 
Long-lived  Vortices  -  K.  Theilhaber  and  C.  K.  Birdsall, 
University  of  California,  Berkeley,  CA  94720 

We  consider  the  equilibrium  and  time-dependent 
behavior  of  the  plasma-wall  sheath  which  is  formed 
when  the  external  magnetic  field  is  parallel  to  the 
wall.  The  presence  of  the  material  boundary  initially 
creates  a  layer  of  charge  depletion  near  the  surface, 
of  thickness  comparable  to  the  ion  gyro-radius .  This 
results  in  a  strongly  nonuniform  electric  field  which 
induces  a  sheared  motion  of  ions  and  electrons 
parallel  to  the  wall .  Physical  arguments  as  well  as 
analytical  and  numerical  studies  suggest  that  this 
flow  is  subject  to  the  Kelvin-Helmholtz  instability, 
with  a  long-term  evolution  Into  a  set  of  vortices 
separating  the  boundary  from  the  plasma.  We  present 
an  analytical  model  for  the  instability,  and  the 
results  of  two-dimensional,  electrostatic  particle 
simulations  of  this  configuration,  in  the  form  of  a 
movie,  showing  initialization  of  the  subsequent 
shear  motion,  and  the  buildup  near  the  wall  of  vortices 
which  live  (at  least)  some  10 's  of  ion  gyroperiods . 


Work  supported  by  DOE,  ONR  and  Varian  Associates/ 
MICRO . 


Contour  plot  of  the  electrostatic 
potential  0  (x,y),  at  U)  t-  100, 
in  a  bounded  plasma  simulation,  using 
the  code  "ES2".  The  confining  walls 
are  ar  x»0  and  x*G4,  the  external 
magnetic  field  is  perpendicular  to 
the  plane  of  the  simulation  region. 

left-hand  potential  depression 
is  a  vortex,  streaming  down  the  wall 
with  a  velocity  of  about  0.1  v  , 
and  the  right-hand  vortex  is  streaming 
up  the  wall , 
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Potential  Drop  and  Transport  in  a  Bounded  Plasma  Including 
Secondary  Electron  Emission  and  Ion  Reflection  at  the  Collector^ 

-  L  A.  SCHWAGERand  C.  K.  BIRDSALL,  Plasma  Theory  and  Stmxtlaiion 
Group,  Electronics  Research  Lab,  University  of  California,  Berkeley  CA 

94720. 

We  simulate  the  model  shown  below  numerically  with  electrostatic  par¬ 
ticle  simulation  using  PDWl  and  analytically  with  a  plasma-sheath  equa¬ 
tion.  Our  numerical  and  ansdytical  results  for  the  surface  potential  and 
plasma  transport  differ  considerably  from  those  using  traditional  Bohm 
sheath  analysis.  For  the  purely  absorbing  surface,  our  analytic  technique 
extends  the  analysis  of  Kuhn^  with  a  temperature  ratio  T,/T,  =  r  =  1  to 
include  an  arbitrary  value  for  r 

Simulation  results  after  equilibration  and  our  time-independent  theory 
agree  exactly  for  r  =  0.1,  1.0,  and  10.0.  With  PDWl  we  observe  that 
initially  drops  from  0  to  4-8  times  the  final  equilibrated  The  minimum 
occurs  at  2-5  times  the  transit  time  of  an  ion  with  the  mean  velocity  after 
the  system  is  initialized.  Note  that  the  presheath  region,  which  we  define 
below  by  d>n,  has  a  zero  electric  field;  hence  the  ‘^resheath  acceleration” 
of  the  ions  occurs  entirely  in  the  source  sheath  within  a  few  Debye  lengths 
lengths  from  the  plasma  source.  Our  potential  results  for  a  large  range  of 
mass  and  temperature  ratios  compare  well  with  Emmert^  who  analytically 
models  the  plasma-sheath  equation  for  a  bounded  system  with  Boltzmann 
electrons. 

When  compared  with  works  such  as  that  of  Chodura*  and  Stangeby^ 
who  rely  on  a  Bohm  type  of  semi-infinite  sheath  model,  out  derivation 
predicts  a  lower  for  a  mass  ratio  of  1836  and  r  =  1;  this  difference  is 
within  20%.  Our  d>w  found  with  r  =  1  and  the  <bu,  derived  by  Stangeby  and 
Chodura  for  r  =  0  agree  to  within  2%.  For  r  >  1  our  results  for  the  ion 
drift  velocity  vd  entering  the  collector  sheath  are  within  4%  of  Stangeby ’s 
assumption  that  this  velocity,. normalized  by  the  ion  thermal  velocity  Vt  is 
(1+ For  T  C  1  our  results  show  amuch  smaller  vd',  e.g.  with  r  »  0.1 
and  a  mass  ratio  of  1836,  vp/vt  —  1.8. 

For  electron-induced  secondary  electron  emission,  we  have  simulated 
emission  beyond  field  reversal  at  the  plate  to  emission  coefficients  F  beyond 
1.  For  emission  values  up  to  the  critical  coefficient  at  field  reversal,  our 
dependence  of  4w  on  F  for  r  =  1  compare  well  to  that  of  Hobbs  and 
Wesson*  who  assume  r  <  1 . 

For  the  reflection  of  plasma  ions,  we  have  simulated  reflection  coeffi¬ 
cients  R  up  to  0.8.  The  wall  potential  decreases  as  we  increase  R  Despite 
an  ion-ion  two-streaming  interaction,  which  generates  substantial  potential 
fluctuations,  the  time-averaged  values  for  and  transport  agree  well  with 
our  time-independent  analysis  which  assumes  r  C  1. 
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